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Abstract 

Using a data sample with a total integrated luminosity of 10.0 pb^^ collected 
at center-of-mass energies of 2.6, 3.07 and 3.65 GeV with BESII, cross sections 
for e~^e~ annihilation into hadronic final states {R values) are measured with 
statistical errors that are smaller than 1%, and systematic errors that are about 
3.5%. The running strong interaction coupling constants al^^(s) and ai^^(M^) 
are determined from the R values. 



1 Introduction 

The R ratio is defined as the lowest level hadronic cross section normalized by the theoretical 
production cross section in e~^e~ annihilation 

p ^ ^ 7* ^ hadrons) 

and it is an important input parameter for precision tests of the Standard Model (SM). The 
errors on R value measurements below 5 GeV have a strong influence on the uncertainties 
of the calculated QED running electromagnetic coupling constant a{s), the muon anomalous 
magnetic moment {g — 2) and global SM fits for the Higgs mass [UlllE]- In addition, precision 
measurements of R values between 2.0 and 3.7 GeV provide a test of perturbative QCD and 
QCD sum rule calculations [ll[5l[6]. 

In 1998 and 1999, R value measurements were made at 91 energy points [TJ [8] between 2 
and 5 GeV by the BESII [9j experiment. The R values were determined from the expression 



R 



where Nf^^^ is the number of observed hadronic events, Ni,g is the number of QED background 
events surviving hadron selection (e+e", z^"*"/^", r+r", 77, etc.), L is the integrated luminosity, 
Etrg is the trigger efficiency for hadronic events, eJlad hadronic efficiency without the 
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simulation of initial state radiation (ISR), and (1 + Sobs) is the effective ISR correction which 
includes the effect of radiated photons on the hadronic acceptance. The average statistical 
errors of the R values are 2 — 4%, and the systematical errors are 5 — 8% depending on the 
energy point; for the latter, the errors associated with the event selection, hadronic efficiency, 
and luminosity are dominant. 

In 2004, large-statistics data samples were accumulated at center-of-mass energies of 2.60, 
3.07 and 3.65 GeV; the total integrated luminosity was 10.0 pb"^ An additional 65.2 nb"^ 
was accumulated at 2.2 GeV for the purpose of tuning the parameters of the hadronic event 
generator. Improvements in the event selection, tuning of generator parameters and luminosity 
measurement have been made in order to decrease the systematic errors. The previously used 
EGS-based (Electron Gamma Shower) p!Q| detector simulation (where hadronic interactions 
were parametrized but not simulated) was replaced by a GEANT3-based one |1H [T2] , where 
both electromagnetic and hadronic interactions are simulated. The consistency between data 
and Monte Carlo (MC) has been validated using many high purity physics channels [13]. With 
these improvements, the errors on the new measured R values are reduced to about 3.5%. 

In this letter, improvements compared to the previous measurements are described, and 
the error analysis and results are reported. Finally, the strong running coupling constant 

(s) and ai^^(M|) are determined from the R values. 

2 Data analysis 

The analysis used for this work is similar to that used in the previous BESII R measurements 
[71 [H] . Two large sources of error in the measurement arise from the event selection and the 
determination of the hadronic efficiency; these are strongly correlated. 

2.1 Selection of hadronic events 

In the BEPC energy region, data collected include processes that originate from beam-beam 
collisions, e~^e~ — > e"*"e^, fi^fi^, t^t^, 77, e~^e~X {X means any possible final state), and 
hadrons (including continuum and resonant states), as well as beam associated backgrounds. 
The observed final state charged particles are e, n, vr, K, and p. To test the hadron selection 
criteria, described below, different types of backgrounds are selected using specialized criteria, 
and most of them are rejected with good efficiency (TJ \S[ I14j. 

The candidate hadronic events are classified by their number of charged tracks. The 
selection of hadronic events is done in two successive steps: one at the track level and the 
other at the event level. 

(I) Track level: only charged tracks that are well fitted to a helix are considered; the point 
of closest approach of the track, signified by {Vx, Vy,Vz), must be within 2 cm of the beam- line 
in the x-y plane with no restriction on Vz ; the angle between the track and the z axis is limited 
by the coverage of the main drift chamber (MDC) to be within | cos^| < 0.84; a momentum 

cut p < Eh{l + 0.1^1 + E"^) removes tracks with unphysically high momentum (Eh is the 
beam energy in units of GeV); a time of fiight (TOF) requirement tTOF < tp + 2 ns is applied 
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to reject events with unphysical times {tp is the expected value for a proton with momentum 
p); the energy deposited in the barrel shower counter (BSC), Esse, must be less than the 
minimum of 1 GeV and O.GEb', and the number of hit layers in the muon counter must be 
smaller than 3. 

A neutral cluster in the BSC is considered to be a photon candidate if the angle between 
the nearest charged track and the cluster is greater than 25°; the difference between the angle 
of the cluster development direction and the photon emission direction in the BSC is less than 
30°; and the number of hit layers in the BSC is larger than two. 

(II) Event level: the requirement that the total deposited energy, E^^q, is greater than 
the maximum of 0.5 GeV and 0.28-Eb eliminates most of the beam-associated backgrounds. 
The selected tracks must not all point into the forward (cos > 0) or the backward (cos 9 < 0) 
hemisphere. For multi-track events {rich ^ 3), no further requirement is used. For two-track 
events, the two charged tracks cannot be back-to-back, so the angle between them is required 
to be less than 165°, and there must be at least two isolated photons with energy E^ > 0.1 
GeV that are well separated from charged tracks, i.e. the distance between the neutral and 
charged tracks at the first layer of the BSC must be larger than 34 cm in the x-y plane, or 
larger than 60 cm in the z direction. 

The above described selection criteria for hadronic events with rich ^ 2 are almost the 
same as used in the previous measurements [3 IB]- In the BEPC energy region, the number of 
events with one observed/reconstructed charged track in BESII accounts for about 8-13% of all 
hadronic events. The omission of one-track (and zero-track) hadronic events introduces some 
uncertainty in the tuning the parameters of the hadronic event generator parameters; this in 
turn induces a sizable systematic error in the hadronic efficiency. However, for single-track 
events, contamination from beam-associated backgrounds is significant. Therefore, a more 
restrictive hadronic event selection is applied |14j: the event must have one charged track that 
is well fitted to a helix (the event can have additional any number of charged tracks with poor 
helix fits). If this charged track is identified as e^, the event is rejected. If its momentum is 
in the range p > 1 GeV and the number of hit layers in the muon counter is larger than 1, the 
event is also rejected. For each single-track event, the number of photons with energy E^ > 0.1 
GeV should be two or more. To further suppress background from fake photons, only events 
with one well fitted charged track and at least one reconstructed vr'' are considered as single- 
track hadronic events. A 1-C fit is applied under the vr'^ — > 77 hypothesis. For candidates 
with more photons (n^ > 3), the 77 pair combination with the smallest is chosen. The 
probability for the 1-C fit is required to be larger than 1%. Figure EJh) shows the invariant 
77 mass distribution for one-track events satisfying requirements for data and MC (the MC is 
normalized to the integrated luminosity of the data) . The good agreement indicates that the 
rate of one-track events in data and MC are commensurate, and the residual beam-associated 
background is small. 

For selected events, the weighted average vertex position, Vz, is determined using 
where Vz{i) is the z coordinate of the ith selected track, is the track fitting chi-square, 
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Xi = Xi /''T'D.o.F., no.o.F. is the number of degrees of freedom, and n is the number of tracks 
in the event. Figure [T] shows the Vz distributions for candidate hadronic events (including 
the residual beam-associated and QED backgrounds) at 2.6, 3.07 and 3.65 GeV. Signal events 
produced by e~^e~ collisions originate near the collision point (in the neighborhood of z = 0), 
and the non-beam-beam backgrounds, such as those from beam-gas and beam-wall scattering, 
are distributed all along the beam direction (In order to show the shape of the Vz distribution 
of residual beam-associated backgrounds clearly, a logarithmic vertical scale is used). The 
number of observed hadronic events iV^^^ is determined by fitting the Vz distribution with 
a Gaussian to describe the hadronic events and an m-th degree polynomial for the residual 
beam associated backgrounds. 

The numbers of residual QED background events, N^g in Eq.([2]), are determined from MC 
simulations using QED event generators with an accuracy of 1% [15j, where 

^bg = L[eee(Tee + ^tJ.^J.<^^J,^l + ^tt<^tt + e^^a^^]. (4) 

Here o"ee is the production cross section for Bhabha events given by the corresponding genera- 
tor, €ee is the efficiency for Bhabha events that pass the hadronic event selection criteria, and 
other symbols have corresponding meanings. The values of €ee and e^^ are about 5 x 10~^, and 
err is 36.45% at 3.65 GeV. The errors on N^g are given in Table El The amount of background 
from e^e~ e~^e~X that survives hadron selection is much smaller than 1% of Nf^g and is 
neglected. 

2.2 Tuning the LUARLW parameters 

The hadronic efficiency is determined using the LUARLW hadronic event generator [16j. 
The physical basis of LUARLW is the Lund area law [17J. The production of hadrons is 
described as the fragmentation of a semi-classical relativistic string [18], and the quark com- 
ponents of the string and decays of unstable particles are handled by subroutines in JESTSET 
[191 [201 El]. 

Both generators LUARLW and JETSET have some phenomenological parameters that 
have to be determined from data. The important parameters in JETSET are PARJ(1 — 
3), which are responsible for the suppression of diquark-diantiquark pair production, quark 
pair production, and the extra suppression of strange diquark production respectively, and 
PARJ{11 — 17), which control the relative probabilities of different spin mesons. The values 
of these parameters have been tuned with LEP data, as well as with information from lower 
energy e~^e~ colliders [22j . The generator LUARLW has also been tuned with BES data taken 
at 2.2, 2.6, 3.07 and 3.65 GeV. The tuned parameters for the Lund area law are mainly the 
dynamical parameter b and those related to the initial multiplicity distributions (including 
neutral clusters and charged tracks) from the string decay p!6] • 

The basic method is to find a set of parameters that make various distributions (especially 
those related to the hadronic selection) simulated by MC agree well with experimental data 
at all of the measured energy points. The distributions used for the data— MC comparison 
are: the multiplicities of charged tracks and neutral clusters, the Vx-y and Vz coordinates of 
charged tracks, the charged track momentum, the polar-angle 9 between tracks and the beam 
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direction, the deposited energy in the BSC, the time of flight, and fractions of vr^, K^, and 
some other short-Hved particles (vr'', Ks, (p, A), etc.. With these distributions, the systematic 
errors corresponding to each criteria used in the hadronic event selection can be determined. 
Figure E] shows, for example, some comparisons between data and the LUARLW MC at 3.65 
GeV, where reasonable agreement is evident. More distributions at other energy points are 
also compared, which can be found in Refs. [THUS]. 

The parameters for detector simulation are studied using data. The constant files for dead 
and hot channels, and the detailed detector responses are inputs for detector simulations. 

2.3 Trigger, luminosity and ISR 

The trigger conditions are almost the same as those used for the R measurements reported in 
Refs. [71 [8], and the details about the values and errors of the trigger efficiency are described 
in Refs. [21]. Since one-track hadronic events are also included in this measurement, the 
TOF back-to-back hit trigger requirement is not used, thereby making the trigger conditions 
somewhat loser than before. The trigger table used in data taking is given in Ref. [Mj. The 
trigger efficiencies etrg for hadronic events is determined as 99.8%, and its associated errors is 
conservatively estimated to be 0.5%. 

The integrated luminosity L is measured with wide-angle (| cos^| < 0.6) Bhabha events. 
The measurement method is similar to that described in Refs. [3 |8]. The Bhabha events are 
selected using only BSC information, and the simulation of the BSC is significantly improved 
by the package (SIMBES) described in [13], which provides better consistency between MC and 
data. In addition, the BSC selection efficiencies, determined by simulation, are corrected using 
correction factors determined with another Bhabha sample, that is selected using only MDC 
information. The efficiency correction factors range from 0.994 to 1.026 with an uncertainty 
of about 1.4% for the different energy points. In addition, the contribution from e~^e~ — > 77 
is subtracted explicitly. As a result, the precision of the luminosity determined is significantly 
improved, and the systematic uncertainties are about 2%. For example, at 2.6 GeV the 
systematic error is 1.9%, of which the trigger efficiency contributes 0.5%, the MC generator 
0.5%, event selection 1.2%, and the Bhabha correction 1.3%. 

An 0{a^) Feynman-diagram-based calculation for the initial state radiative (ISR) cor- 
rection is used in both the calculation of the ISR factor (1 + 6obs) and the simulation of 
radiative events by LUARLW. A detailed description of the ISR treatment can be found in 
Refs. [25l [26l [271 [28l [29] . the ISR simulations and calculations, the contributions from both 
continuum and resonances are considered. The inclusive hadronic cross section below 5 GeV 
uses the experimental value, and above 5 GeV uses the theoretical value predicted by QCD 
|30j . The quantities related to the narrow J/'i/; and ■0' are treated analytically. 

Figure E] shows the detection efficiencies for hadronic events simulated with LUARLW for 
the cases where an initial-state e"*" or e~ radiates a photon with energy fraction k = E^/Ef, 
and those events with n^h > 1 are selected. Curves fitted to e{k) are used in the calculation 
of (1 + Sobs) at different energies. The values of (1 + Sobs) and their errors are listed in Tables 
[Hand [21 respectively. 
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3 Error analysis 



The Feynman-diagram-based ISR simulated angle and momentum distributions for the 
radiated photon are built into LUARLW, which allows the hadronic efficiency Chad including 
radiative effects to be obtained, averaging over the ISR spectrum. The number of hadronic 
events A'^^^d' hadronic efficiency ehad, and their errors are correlated. The equivalent 
number of hadronic events, which corresponds to the number of hadronic events produced at 
the collision point, is defined as 

ATobs 

Nhad = (5) 
^had 

The combined systematic error associated with the event selection and hadronic efficiency is 
denoted as ANhad- This error is caused by discrepancies between data and MC samples for 
the hadronic selection criteria discussed in Section 2.1. In addition, an uncertainty associated 
with the parameters of the MC hadronization model is estimated to be about 1% by comparing 
different sets of tuned parameters, and is included in the error of the hadronic event efficiency. 

The error on A^^^^ due to the choice of the degree of the polynomial used in the fitting 
is less than 0.7%. The fit errors for Nf^J^, which are calculated from the uncertainties in the 
fitted parameters of the Gaussian signal peaks, are 1.34% at 2.6 GeV, 1.11% at 3.07 GeV, and 
0.73% at 3.65 GeV. The total AN^ad is the quadratic sum of all fractional errors. 

The uncertainty in the effective ISR factor (1 + 6obs) due to errors of the hadronic cross 
sections at the different effective energies for radiative events is considered (the errors on the 
hadronic cross section given in the PDG06 tables [30] are used); these decrease with increasing 
energy from 0.9% to 0.1%. For comparison, another approach based on structure functions 
[31j is also used at all energy points. The differences in the values of (1 + 5) for the two schemes 
are smaller than 1.1%. 

A conclusion of the KLN theorem is that the radiative corrections due to final state radia- 
tion (FSR) are negligible for a measurement of the inclusive hadronic cross section that sums 
over all hadronic final states [32]. At the present level of precision, the FSR correction factor 
in Eq.([2]) can be neglected. However, the absence of final state radiation in the event generator 
introduces some error into the determination of the hadronic event detection efficiency. The 
masses of the produced hadrons in the final states are much greater than that of the initial 
radiative e^. As a result, the effect of FSR is much weaker than initial bremsstrahlung. Its 
influence is estimated to be 0.5% and is included in the error. 

The 0-track hadronic events are not selected in this analysis, and the influence of 0-track 
events on the parameter tuning of LUARLW is not considered. This introduces some error 
into the hadronic event efficiency determination. Events with no charged tracks cannot be well 
separated from background. The fraction of 0-track events is estimated from the MC to be 
3.4% at 2.6 GeV, 2.9% at 3.07 GeV, and 2.4% at 3.65 GeV. If the difference for 0-track events 
between MC and data is conservatively assumed to be 20%, the errors for the lost /unobserved 
0-track events are 0.7%, 0.6% and 0.5%, respectively. The error related to 0-track events is 
included into the error of N^ad deflned in Eq. ([5]). 

In this analysis, hadronic events are classifled according to their number of charged tracks. 
Therefore, errors in the tracking efficiency Cj^fc) the differences in the track reconstruction 
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efficiency between data and MC, introduce some error into the classification and counting of 
the number of events. For an event with rich charged tracks, the probability that of rich 
tracks are wrongly constructed roughly obeys a binomial distribution B{ncr]nch^(^trk)^ where 
the parameter atrk ~ 2% is the tracking efficiency error at BESII. The R value measurement is, 
in fact, a counting of the number of hadronic events, so only those cases where all rich tracks 
in an event are incorrectly reconstructed {rier = rich) will introduce an error into ANhad- 
Considering the distribution of charged multiplicity P{nch) for the inclusive hadronic sample 
(such as shown in Figure [2]^ a)), the effective error of tracking efficiency is 

Aetrfc = ^ P (rich) B (rier = rich] rich, CTtrk)- (6) 

The estimated values of Aetrk are listed in Table [2j Since the fraction of single-track events 
decreases with increasing center-of-mass energy, the error Aetrk also decreases with energy. 

The errors on the R values are determined with Eq. ([2]) (including all errors analyzed in 
last section and and estimated by Eq. ([6])). The final R values are 2.18 it 0.02 it 0.08 at 2.6 
GeV, 2. 13 ±0.02 ±0.07 at 3.07 GeV, and 2. 14 ±0.01 ±0.07 at 3.65 GeV, which are summarized 
in Table [H 

As a cross check, the R values are determined using the relation 

'''''''' ~ a^,,Letrg-ehad{l + Sy ^'^ 

where Ihad is the hadronic efficiency averaged over all of the ISR spectrum, and (1 ± 5) is 
the corresponding theoretical ISR factor. The ISR scheme used to simulate Ihad-, e^ad *o 
calculate (1±5) and {l + dobs) are the same in order to keep the consistency between theoretical 
calculation and simulation. The R values determined with Eq. ([7]) are 2.17 ±0.01 ±0.07 at 2.6 
GeV, 2.13 ± 0.01 ± 0.07 at 3.07 GeV, and 2.16 ± 0.01 ± 0.08 at 3.65 GeV. The mean R values 
obtained using Eqs. ([2]) and ([7]) are consistent to within 1%. 

A cross check is also made by selecting hadronic events with rich > 2 as was done in 
Refs. O [8]. In this case, the R values at the three energy points are 2.20 ± 0.02 ± 0.08, 
2.13 ± 0.02 ± 0.07, and 2.15 ± 0.01 ± 0.08, respectively. The differences in the mean R values 
determined by selecting hadronic events with rich ^ 1 and rich ^ 2 are consistent within 1%. 



4 Results and discussion 

Tables [U and [2] list the quantities used in the determination of R using Eq. ^ and 
the contributions to the total error. The results are shown in Fig. HI together with previous 
measurements. The errors on the R values reported here are about 3.5%. The R values are 
consistent within errors with the prediction of perturbative QCD [4j. 

Compared with our previous results [71[8], the measurement precision has been improved 
due to three main refinements to the analysis: (1) the simulation of BES including both of the 
hadronic and electromagnetic interactions with a GEANT3 based package SIMBES that has 
a more detailed geometrical description and matter definition for the sub-detectors; (2) large 
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data samples are taken at each energy point, with statistical errors smaller than 1%; (3) the 
selected hadronic event sample is expanded to include one-track events, which supplies more 
information to the tuning of LUARLW, and results in the improved values of parameters and 
hadronic efficiency. 

In another BESII work, parts of the data sample taken at 3.65 GeV with a luminosity of 
5.536 pb~^ and at 3.665 GeV with a luminosity of 998.2 nb~^ are used, the hadronic events with 
more than 2-tracks {rich ^ 3) are selected, and the averaged R value is R = 2.218ib0.019ib0.089 
which has an error of 4.1% [33j. 

Based on the R values in this work and the perturbative QCD expansion that computes 
Rqcd{o-s) to 0{a^g) [31], the strong running coupling constant ai^^(s) can be determined at 
each energy point [351 136]. The obtained ah (s) values are evolved to 5 GeV, and the weighted 
average of the measurements Q!i'^''(25 GeV^) is listed in Table [3l When evaluated at the Mz 
scale, the resulting value is ai^^(M|) = 0.117j]o;Qj;y, which agrees with the world average value 
within the quoted errors [1]. 



Table 1: Items used in the determination of R at each energy point. 





L (pb-^) 


ATobs 




(%) 


(1 + Sobs) 


R 




^sys 


2.60 


1.222 


24026 


193 


63.81 


1.08 


2.18 


0.02 


0.08 


3.07 


2.291 


33933 


208 


67.63 


1.11 


2.13 


0.02 


0.07 


3.65 


6.485 


83767 


4937 


71.83 


1.21 


2.14 


0.01 


0.07 



Table 2: Summary of the systematic errors (%). 



^cm(GeV) 


L 


Nhad 


Nbg 




^trg 


(1 + Sobs) 


Total 


2.60 


2.00 


2.79 


0.05 


0.32 


0.50 


1.18 


3.68 


3.07 


1.96 


2.53 


0.05 


0.29 


0.50 


1.15 


3.45 


3.65 


1.38 


2.74 


0.35 


0.26 


0.50 


1.10 


3.33 



Table 3: as{s) determined from R values at 2.600, 3.070, and 3.650 GeV, and evolved to 
5 GeV. The first and second errors are statistical and systematic, respectively. Shown 
in the last two columns are the weighted averages of the three measurements at 5 GeV 
and Mz. 



^s{GeY) ^f\^) ^W(25CkV^) aW(25GeV') W^iM^ 



9 en n 9fifi+U.03U+U.125 n oi 9+U.UI8+U.O68 

Z.UU U.ZUU_g 030-0.116 U.ZiZ_o 019-O.O86 

q 07 n 1 Q9+0. 029+0.103 n 1 «Q+0.022+0.074 n. 9nQ+0.044 n 1 1 7+O.OI2 

O.Ui U.lJZ_o. 029-0.101 U.iUy_o.023-0.086 U.ZUJ_o.050 U.ii/_o.017 

Q Rc; n On7+0-015+0-104 n 1 SQ+0-012+0.082 

O.UO U.ZU/_o.015-0.104 U. iOy_o.oi3-0.091 
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Figure 1: Distributions of Vz for candidate hadronic events. The distributions are 
fitted by a Gaussian to represent the signal and a polynomial to describe the residual 
beam-associated backgrounds. 
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Figure 2: The normalized distributions for data (dots with error bars ) and the LUARLW 
MC (histograms) at 3.65 GeV with full detector simulation: (a) multiplicity of charged 
tracks; (b) multiplicity of neutral clusters; (c) momenta p of charged tracks; (d) polar- 
angles between charged tracks and beam direction, cos 6'; (e) deposited energies in the 
BSC; invariant masses of (f) Ks — ^ tt+tt"; (g) A — > p7r~ and (h) 77 decays 

respectively. 
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Figure 3: Radiative efficiencies e{k) for selected hadronic events with rich > 1 at some 
values of k. 
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Figure 4: R values reported here together with prevous measurements below 5 GeV. 
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Abstract 

Using a data sample with a total integrated luminosity of 10.0 pb^^ collected 
at center-of-mass energies of 2.6, 3.07 and 3.65 GeV with BESII, cross sections 
for e~^e~ annihilation into hadronic final states {R values) are measured with 
statistical errors that are smaller than 1%, and systematic errors that are about 
3.5%. The running strong interaction coupling constants al^^(s) and ai^^(M^) 
are determined from the R values. 



1 Introduction 

The R ratio is defined as the lowest level hadronic cross section normalized by the theoretical 
production cross section in e~^e~ annihilation 

p ^ ^ 7* ^ hadrons) 

and it is an important input parameter for precision tests of the Standard Model (SM). The 
errors on R value measurements below 5 GeV have a strong influence on the uncertainties 
of the calculated QED running electromagnetic coupling constant a{s), the muon anomalous 
magnetic moment {g — 2) and global SM fits for the Higgs mass [UlllE]- In addition, precision 
measurements of R values between 2.0 and 3.7 GeV provide a test of perturbative QCD and 
QCD sum rule calculations [ll[5l[6]. 

In 1998 and 1999, R value measurements were made at 91 energy points [TJ [8] between 2 
and 5 GeV by the BESII [9j experiment. The R values were determined from the expression 



R 



where Nf^^^ is the number of observed hadronic events, Ni,g is the number of QED background 
events surviving hadron selection (e+e", z^"*"/^", r+r", 77, etc.), L is the integrated luminosity, 
Etrg is the trigger efficiency for hadronic events, eJlad hadronic efficiency without the 
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simulation of initial state radiation (ISR), and (1 + Sobs) is the effective ISR correction which 
includes the effect of radiated photons on the hadronic acceptance. The average statistical 
errors of the R values are 2 — 4%, and the systematical errors are 5 — 8% depending on the 
energy point; for the latter, the errors associated with the event selection, hadronic efficiency, 
and luminosity are dominant. 

In 2004, large-statistics data samples were accumulated at center-of-mass energies of 2.60, 
3.07 and 3.65 GeV; the total integrated luminosity was 10.0 pb"^ An additional 65.2 nb"^ 
was accumulated at 2.2 GeV for the purpose of tuning the parameters of the hadronic event 
generator. Improvements in the event selection, tuning of generator parameters and luminosity 
measurement have been made in order to decrease the systematic errors. The previously used 
EGS-based (Electron Gamma Shower) p!Q| detector simulation (where hadronic interactions 
were parametrized but not simulated) was replaced by a GEANT3-based one |1H [T2] , where 
both electromagnetic and hadronic interactions are simulated. The consistency between data 
and Monte Carlo (MC) has been validated using many high purity physics channels [13]. With 
these improvements, the errors on the new measured R values are reduced to about 3.5%. 

In this letter, improvements compared to the previous measurements are described, and 
the error analysis and results are reported. Finally, the strong running coupling constant 

(s) and ai^^(M|) are determined from the R values. 

2 Data analysis 

The analysis used for this work is similar to that used in the previous BESII R measurements 
[71 [H] . Two large sources of error in the measurement arise from the event selection and the 
determination of the hadronic efficiency; these are strongly correlated. 

2.1 Selection of hadronic events 

In the BEPC energy region, data collected include processes that originate from beam-beam 
collisions, e~^e~ — > e"*"e^, fi^fi^, t^t^, 77, e~^e~X {X means any possible final state), and 
hadrons (including continuum and resonant states), as well as beam associated backgrounds. 
The observed final state charged particles are e, n, vr, K, and p. To test the hadron selection 
criteria, described below, different types of backgrounds are selected using specialized criteria, 
and most of them are rejected with good efficiency (TJ \S[ I14j. 

The candidate hadronic events are classified by their number of charged tracks. The 
selection of hadronic events is done in two successive steps: one at the track level and the 
other at the event level. 

(I) Track level: only charged tracks that are well fitted to a helix are considered; the point 
of closest approach of the track, signified by {Vx, Vy,Vz), must be within 2 cm of the beam- line 
in the x-y plane with no restriction on Vz ; the angle between the track and the z axis is limited 
by the coverage of the main drift chamber (MDC) to be within | cos^| < 0.84; a momentum 

cut p < Eh{l + 0.1^1 + E"^) removes tracks with unphysically high momentum (Eh is the 
beam energy in units of GeV); a time of fiight (TOF) requirement tTOF < tp + 2 ns is applied 
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to reject events with unphysical times {tp is the expected value for a proton with momentum 
p); the energy deposited in the barrel shower counter (BSC), Esse, must be less than the 
minimum of 1 GeV and O.GEb', and the number of hit layers in the muon counter must be 
smaller than 3. 

A neutral cluster in the BSC is considered to be a photon candidate if the angle between 
the nearest charged track and the cluster is greater than 25°; the difference between the angle 
of the cluster development direction and the photon emission direction in the BSC is less than 
30°; and the number of hit layers in the BSC is larger than two. 

(II) Event level: the requirement that the total deposited energy, E^^q, is greater than 
the maximum of 0.5 GeV and 0.28-Eb eliminates most of the beam-associated backgrounds. 
The selected tracks must not all point into the forward (cos > 0) or the backward (cos 9 < 0) 
hemisphere. For multi-track events {rich ^ 3), no further requirement is used. For two-track 
events, the two charged tracks cannot be back-to-back, so the angle between them is required 
to be less than 165°, and there must be at least two isolated photons with energy E^ > 0.1 
GeV that are well separated from charged tracks, i.e. the distance between the neutral and 
charged tracks at the first layer of the BSC must be larger than 34 cm in the x-y plane, or 
larger than 60 cm in the z direction. 

The above described selection criteria for hadronic events with rich ^ 2 are almost the 
same as used in the previous measurements [3 IB]- In the BEPC energy region, the number of 
events with one observed/reconstructed charged track in BESII accounts for about 8-13% of all 
hadronic events. The omission of one-track (and zero-track) hadronic events introduces some 
uncertainty in the tuning the parameters of the hadronic event generator parameters; this in 
turn induces a sizable systematic error in the hadronic efficiency. However, for single-track 
events, contamination from beam-associated backgrounds is significant. Therefore, a more 
restrictive hadronic event selection is applied |14j: the event must have one charged track that 
is well fitted to a helix (the event can have additional any number of charged tracks with poor 
helix fits). If this charged track is identified as e^, the event is rejected. If its momentum is 
in the range p > 1 GeV and the number of hit layers in the muon counter is larger than 1, the 
event is also rejected. For each single-track event, the number of photons with energy E^ > 0.1 
GeV should be two or more. To further suppress background from fake photons, only events 
with one well fitted charged track and at least one reconstructed vr'' are considered as single- 
track hadronic events. A 1-C fit is applied under the vr'^ — > 77 hypothesis. For candidates 
with more photons (n^ > 3), the 77 pair combination with the smallest is chosen. The 
probability for the 1-C fit is required to be larger than 1%. Figure EJh) shows the invariant 
77 mass distribution for one-track events satisfying requirements for data and MC (the MC is 
normalized to the integrated luminosity of the data) . The good agreement indicates that the 
rate of one-track events in data and MC are commensurate, and the residual beam-associated 
background is small. 

For selected events, the weighted average vertex position, Vz, is determined using 
where Vz{i) is the z coordinate of the ith selected track, is the track fitting chi-square, 
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Xi = Xi /''T'D.o.F., no.o.F. is the number of degrees of freedom, and n is the number of tracks 
in the event. Figure [T] shows the Vz distributions for candidate hadronic events (including 
the residual beam-associated and QED backgrounds) at 2.6, 3.07 and 3.65 GeV. Signal events 
produced by e~^e~ collisions originate near the collision point (in the neighborhood of z = 0), 
and the non-beam-beam backgrounds, such as those from beam-gas and beam-wall scattering, 
are distributed all along the beam direction (In order to show the shape of the Vz distribution 
of residual beam-associated backgrounds clearly, a logarithmic vertical scale is used). The 
number of observed hadronic events iV^^^ is determined by fitting the Vz distribution with 
a Gaussian to describe the hadronic events and an m-th degree polynomial for the residual 
beam associated backgrounds. 

The numbers of residual QED background events, N^g in Eq.([2]), are determined from MC 
simulations using QED event generators with an accuracy of 1% [15j, where 

^bg = L[eee(Tee + ^tJ.^J.<^^J,^l + ^tt<^tt + e^^a^^]. (4) 

Here o"ee is the production cross section for Bhabha events given by the corresponding genera- 
tor, €ee is the efficiency for Bhabha events that pass the hadronic event selection criteria, and 
other symbols have corresponding meanings. The values of €ee and e^^ are about 5 x 10~^, and 
err is 36.45% at 3.65 GeV. The errors on N^g are given in Table El The amount of background 
from e^e~ e~^e~X that survives hadron selection is much smaller than 1% of Nf^g and is 
neglected. 

2.2 Tuning the LUARLW parameters 

The hadronic efficiency is determined using the LUARLW hadronic event generator [16j. 
The physical basis of LUARLW is the Lund area law [17J. The production of hadrons is 
described as the fragmentation of a semi-classical relativistic string [18], and the quark com- 
ponents of the string and decays of unstable particles are handled by subroutines in JESTSET 
[191 [201 El]. 

Both generators LUARLW and JETSET have some phenomenological parameters that 
have to be determined from data. The important parameters in JETSET are PARJ(1 — 
3), which are responsible for the suppression of diquark-diantiquark pair production, quark 
pair production, and the extra suppression of strange diquark production respectively, and 
PARJ{11 — 17), which control the relative probabilities of different spin mesons. The values 
of these parameters have been tuned with LEP data, as well as with information from lower 
energy e~^e~ colliders [22j . The generator LUARLW has also been tuned with BES data taken 
at 2.2, 2.6, 3.07 and 3.65 GeV. The tuned parameters for the Lund area law are mainly the 
dynamical parameter b and those related to the initial multiplicity distributions (including 
neutral clusters and charged tracks) from the string decay p!6] • 

The basic method is to find a set of parameters that make various distributions (especially 
those related to the hadronic selection) simulated by MC agree well with experimental data 
at all of the measured energy points. The distributions used for the data— MC comparison 
are: the multiplicities of charged tracks and neutral clusters, the Vx-y and Vz coordinates of 
charged tracks, the charged track momentum, the polar-angle 9 between tracks and the beam 
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direction, the deposited energy in the BSC, the time of flight, and fractions of vr^, K^, and 
some other short-Hved particles (vr'', Ks, (p, A), etc.. With these distributions, the systematic 
errors corresponding to each criteria used in the hadronic event selection can be determined. 
Figure E] shows, for example, some comparisons between data and the LUARLW MC at 3.65 
GeV, where reasonable agreement is evident. More distributions at other energy points are 
also compared, which can be found in Refs. [THUS]. 

The parameters for detector simulation are studied using data. The constant files for dead 
and hot channels, and the detailed detector responses are inputs for detector simulations. 

2.3 Trigger, luminosity and ISR 

The trigger conditions are almost the same as those used for the R measurements reported in 
Refs. [71 [8], and the details about the values and errors of the trigger efficiency are described 
in Refs. [21]. Since one-track hadronic events are also included in this measurement, the 
TOF back-to-back hit trigger requirement is not used, thereby making the trigger conditions 
somewhat loser than before. The trigger table used in data taking is given in Ref. [Mj. The 
trigger efficiencies etrg for hadronic events is determined as 99.8%, and its associated errors is 
conservatively estimated to be 0.5%. 

The integrated luminosity L is measured with wide-angle (| cos^| < 0.6) Bhabha events. 
The measurement method is similar to that described in Refs. [3 |8]. The Bhabha events are 
selected using only BSC information, and the simulation of the BSC is significantly improved 
by the package (SIMBES) described in [13], which provides better consistency between MC and 
data. In addition, the BSC selection efficiencies, determined by simulation, are corrected using 
correction factors determined with another Bhabha sample, that is selected using only MDC 
information. The efficiency correction factors range from 0.994 to 1.026 with an uncertainty 
of about 1.4% for the different energy points. In addition, the contribution from e~^e~ — > 77 
is subtracted explicitly. As a result, the precision of the luminosity determined is significantly 
improved, and the systematic uncertainties are about 2%. For example, at 2.6 GeV the 
systematic error is 1.9%, of which the trigger efficiency contributes 0.5%, the MC generator 
0.5%, event selection 1.2%, and the Bhabha correction 1.3%. 

An 0{a^) Feynman-diagram-based calculation for the initial state radiative (ISR) cor- 
rection is used in both the calculation of the ISR factor (1 + 6obs) and the simulation of 
radiative events by LUARLW. A detailed description of the ISR treatment can be found in 
Refs. [25l [26l [271 [28l [29] . the ISR simulations and calculations, the contributions from both 
continuum and resonances are considered. The inclusive hadronic cross section below 5 GeV 
uses the experimental value, and above 5 GeV uses the theoretical value predicted by QCD 
|30j . The quantities related to the narrow J/'i/; and ■0' are treated analytically. 

Figure E] shows the detection efficiencies for hadronic events simulated with LUARLW for 
the cases where an initial-state e"*" or e~ radiates a photon with energy fraction k = E^/Ef, 
and those events with n^h > 1 are selected. Curves fitted to e{k) are used in the calculation 
of (1 + Sobs) at different energies. The values of (1 + Sobs) and their errors are listed in Tables 
[Hand [21 respectively. 
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3 Error analysis 



The Feynman-diagram-based ISR simulated angle and momentum distributions for the 
radiated photon are built into LUARLW, which allows the hadronic efficiency Chad including 
radiative effects to be obtained, averaging over the ISR spectrum. The number of hadronic 
events A'^^^d' hadronic efficiency ehad, and their errors are correlated. The equivalent 
number of hadronic events, which corresponds to the number of hadronic events produced at 
the collision point, is defined as 

ATobs 

Nhad = (5) 
^had 

The combined systematic error associated with the event selection and hadronic efficiency is 
denoted as ANhad- This error is caused by discrepancies between data and MC samples for 
the hadronic selection criteria discussed in Section 2.1. In addition, an uncertainty associated 
with the parameters of the MC hadronization model is estimated to be about 1% by comparing 
different sets of tuned parameters, and is included in the error of the hadronic event efficiency. 

The error on A^^^^ due to the choice of the degree of the polynomial used in the fitting 
is less than 0.7%. The fit errors for Nf^J^, which are calculated from the uncertainties in the 
fitted parameters of the Gaussian signal peaks, are 1.34% at 2.6 GeV, 1.11% at 3.07 GeV, and 
0.73% at 3.65 GeV. The total AN^ad is the quadratic sum of all fractional errors. 

The uncertainty in the effective ISR factor (1 + 6obs) due to errors of the hadronic cross 
sections at the different effective energies for radiative events is considered (the errors on the 
hadronic cross section given in the PDG06 tables [30] are used); these decrease with increasing 
energy from 0.9% to 0.1%. For comparison, another approach based on structure functions 
[31j is also used at all energy points. The differences in the values of (1 + 5) for the two schemes 
are smaller than 1.1%. 

A conclusion of the KLN theorem is that the radiative corrections due to final state radia- 
tion (FSR) are negligible for a measurement of the inclusive hadronic cross section that sums 
over all hadronic final states [32]. At the present level of precision, the FSR correction factor 
in Eq.([2]) can be neglected. However, the absence of final state radiation in the event generator 
introduces some error into the determination of the hadronic event detection efficiency. The 
masses of the produced hadrons in the final states are much greater than that of the initial 
radiative e^. As a result, the effect of FSR is much weaker than initial bremsstrahlung. Its 
influence is estimated to be 0.5% and is included in the error. 

The 0-track hadronic events are not selected in this analysis, and the influence of 0-track 
events on the parameter tuning of LUARLW is not considered. This introduces some error 
into the hadronic event efficiency determination. Events with no charged tracks cannot be well 
separated from background. The fraction of 0-track events is estimated from the MC to be 
3.4% at 2.6 GeV, 2.9% at 3.07 GeV, and 2.4% at 3.65 GeV. If the difference for 0-track events 
between MC and data is conservatively assumed to be 20%, the errors for the lost /unobserved 
0-track events are 0.7%, 0.6% and 0.5%, respectively. The error related to 0-track events is 
included into the error of N^ad deflned in Eq. ([5]). 

In this analysis, hadronic events are classifled according to their number of charged tracks. 
Therefore, errors in the tracking efficiency Cj^fc) the differences in the track reconstruction 
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efficiency between data and MC, introduce some error into the classification and counting of 
the number of events. For an event with rich charged tracks, the probability that of rich 
tracks are wrongly constructed roughly obeys a binomial distribution B{ncr]nch^(^trk)^ where 
the parameter atrk ~ 2% is the tracking efficiency error at BESII. The R value measurement is, 
in fact, a counting of the number of hadronic events, so only those cases where all rich tracks 
in an event are incorrectly reconstructed {rier = rich) will introduce an error into ANhad- 
Considering the distribution of charged multiplicity P{nch) for the inclusive hadronic sample 
(such as shown in Figure [2]^ a)), the effective error of tracking efficiency is 

Aetrfc = ^ P (rich) B (rier = rich] rich, CTtrk)- (6) 

The estimated values of Aetrk are listed in Table [2j Since the fraction of single-track events 
decreases with increasing center-of-mass energy, the error Aetrk also decreases with energy. 

The errors on the R values are determined with Eq. ([2]) (including all errors analyzed in 
last section and and estimated by Eq. ([6])). The final R values are 2.18 it 0.02 it 0.08 at 2.6 
GeV, 2. 13 ±0.02 ±0.07 at 3.07 GeV, and 2. 14 ±0.01 ±0.07 at 3.65 GeV, which are summarized 
in Table [H 

As a cross check, the R values are determined using the relation 

'''''''' ~ a^,,Letrg-ehad{l + Sy ^'^ 

where Ihad is the hadronic efficiency averaged over all of the ISR spectrum, and (1 ± 5) is 
the corresponding theoretical ISR factor. The ISR scheme used to simulate Ihad-, e^ad *o 
calculate (1±5) and {l + dobs) are the same in order to keep the consistency between theoretical 
calculation and simulation. The R values determined with Eq. ([7]) are 2.17 ±0.01 ±0.07 at 2.6 
GeV, 2.13 ± 0.01 ± 0.07 at 3.07 GeV, and 2.16 ± 0.01 ± 0.08 at 3.65 GeV. The mean R values 
obtained using Eqs. ([2]) and ([7]) are consistent to within 1%. 

A cross check is also made by selecting hadronic events with rich > 2 as was done in 
Refs. O [8]. In this case, the R values at the three energy points are 2.20 ± 0.02 ± 0.08, 
2.13 ± 0.02 ± 0.07, and 2.15 ± 0.01 ± 0.08, respectively. The differences in the mean R values 
determined by selecting hadronic events with rich ^ 1 and rich ^ 2 are consistent within 1%. 



4 Results and discussion 

Tables [U and [2] list the quantities used in the determination of R using Eq. ^ and 
the contributions to the total error. The results are shown in Fig. HI together with previous 
measurements. The errors on the R values reported here are about 3.5%. The R values are 
consistent within errors with the prediction of perturbative QCD [4j. 

Compared with our previous results [71[8], the measurement precision has been improved 
due to three main refinements to the analysis: (1) the simulation of BES including both of the 
hadronic and electromagnetic interactions with a GEANT3 based package SIMBES that has 
a more detailed geometrical description and matter definition for the sub-detectors; (2) large 
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data samples are taken at each energy point, with statistical errors smaller than 1%; (3) the 
selected hadronic event sample is expanded to include one-track events, which supplies more 
information to the tuning of LUARLW, and results in the improved values of parameters and 
hadronic efficiency. 

In another BESII work, parts of the data sample taken at 3.65 GeV with a luminosity of 
5.536 pb~^ and at 3.665 GeV with a luminosity of 998.2 nb~^ are used, the hadronic events with 
more than 2-tracks {rich ^ 3) are selected, and the averaged R value is R = 2.218ib0.019ib0.089 
which has an error of 4.1% [33j. 

Based on the R values in this work and the perturbative QCD expansion that computes 
Rqcd{o-s) to 0{a^g) [31], the strong running coupling constant ai^^(s) can be determined at 
each energy point [351 136]. The obtained ah (s) values are evolved to 5 GeV, and the weighted 
average of the measurements Q!i'^''(25 GeV^) is listed in Table [3l When evaluated at the Mz 
scale, the resulting value is ai^^(M|) = 0.117j]o;Qj;y, which agrees with the world average value 
within the quoted errors [1]. 



Table 1: Items used in the determination of R at each energy point. 





L (pb-^) 


ATobs 




(%) 


(1 + Sobs) 


R 




^sys 


2.60 


1.222 


24026 


193 


63.81 


1.08 


2.18 


0.02 


0.08 


3.07 


2.291 


33933 


208 


67.63 


1.11 


2.13 


0.02 


0.07 


3.65 


6.485 


83767 


4937 


71.83 


1.21 


2.14 


0.01 


0.07 



Table 2: Summary of the systematic errors (%). 



^cm(GeV) 


L 


Nhad 


Nbg 




^trg 


(1 + Sobs) 


Total 


2.60 


2.00 


2.79 


0.05 


0.32 


0.50 


1.18 


3.68 


3.07 


1.96 


2.53 


0.05 


0.29 


0.50 


1.15 


3.45 


3.65 


1.38 


2.74 


0.35 


0.26 


0.50 


1.10 


3.33 



Table 3: as{s) determined from R values at 2.600, 3.070, and 3.650 GeV, and evolved to 
5 GeV. The first and second errors are statistical and systematic, respectively. Shown 
in the last two columns are the weighted averages of the three measurements at 5 GeV 
and Mz. 



^s{GeY) ^f\^) ^W(25CkV^) aW(25GeV') W^iM^ 



9 en n 9fifi+U.03U+U.125 n oi 9+U.UI8+U.O68 

Z.UU U.ZUU_g 030-0.116 U.ZiZ_o 019-O.O86 

q 07 n 1 Q9+0. 029+0.103 n 1 «Q+0.022+0.074 n. 9nQ+0.044 n 1 1 7+O.OI2 

O.Ui U.lJZ_o. 029-0.101 U.iUy_o.023-0.086 U.ZUJ_o.050 U.ii/_o.017 

Q Rc; n On7+0-015+0-104 n 1 SQ+0-012+0.082 

O.UO U.ZU/_o.015-0.104 U. iOy_o.oi3-0.091 
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Figure 1: Distributions of Vz for candidate hadronic events. The distributions are 
fitted by a Gaussian to represent the signal and a polynomial to describe the residual 
beam-associated backgrounds. 
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Figure 2: The normalized distributions for data (dots with error bars ) and the LUARLW 
MC (histograms) at 3.65 GeV with full detector simulation. From left to right and upper 
to down are: (a) multiplicity of charged tracks; (b) multiplicity of neutral clusters; 
(c) momenta p of charged tracks; (d) polar-angles between charged tracks and beam 
direction, cos 6'; (e) deposited energies in the BSC; invariant masses of (f) Ks tt+tt"; 
(g) A piT' and (h) 7r° 77 decays respectively. 



e{k) 



■ Ecm=2.60 GeV 
^ Ecnn=3.07 GeV 
- Ecnn=3.65 GeV 



O I ' ^ ' 1 

O 0.25 0.5 0.75 1 

Figure 3: Radiative efficiencies e{k) for selected hadronic events with rich > 1 at some 
values of k. 
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Figure 4: R values reported here together with prevous measurements below 5 GeV. 
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